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METHOD AND INSTALLATION FOR OPTO -ACOUSTICAL IMAGERY 

The present invention relates to a method and 
installation for acousto-optical imaging. 

In acousto-optical imaging, one detects the beat of a 
local oscillator with an acoustic component of a signal 
wave scattered by an object to be imaged, shifted in 
frequency by the vibration at an acoustic frequency of 
a point of said object to be imaged from which one 
seeks to obtain information of an optical nature. 

In monopixel acousto-optical imaging, the measurement 
of the beat between the component of the signal wave 
scattered without frequency shift, which serves as 
local oscillator, and the acoustic component of the 
signal wave, which carries the information, exhibits 
significant noise since the measurement of said beat is 
performed at just one point of the detection plane 
whereas each of these two components varies in a random 
manner in this plane. It is necessary to perform a 
summation over time of the square of the amplitude of 
said beat to obtain information exhibiting a better 
signal/noise ratio . 

To alleviate this drawback, it is possible to use a 
multipixel detection device, by performing a summation 
over the pixels of the detection device, rather than in 
time, such as described in "Ultrasonic tagging of 
photon paths in scattering media: parallel speckle 
modulation processing" , Optics Letters, Vol. 24, No. 3, 
1st February 1999, page 181. In this context, it is 
necessary to modulate the power of the laser at a 
frequency close to the acoustic frequency, so that the 
beat between the acoustic component of the signal wave, 
carrying the information, and the modulation lateral 
band (the local oscillator) is of sufficiently low 
frequency to be detected by a multipixel detection 
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device, which in general possesses a low acquisition 
frequency. Nevertheless, a major problem remains, in 
that the weight of the local oscillator is in general 
too low. The heterodyne gain is then too low to be able 
5 to perform heterodyne detection with optimal noise. 

The present invention is aimed in particular at 
alleviating these drawbacks. 

10 For this purpose, there is provided according to the 
invention, a method of acousto-optical imaging of an 
object to be imaged comprising the steps consisting in: 

(a) generating an incident optical wave, of 
frequency fi, and a reference optical wave, of frequency 

15 f R , this reference wave being coherent with the incident 
wave, and exhibiting therewith a known phase difference 
0i(t) , 

(b) vibrating in a first object direction and at 
an acoustic frequency f A , a zone of the object to be 

2 0 imaged with the aid of a vibration generating device, 

(c) applying said incident wave to the object to 
be imaged, and thus generating a scattered signal wave, 

(d) applying at least a part of the scattered 
signal wave to a detection device, 

25 (e) applying the reference wave to the detection 

device without making it pass through the object to be 
imaged, thereby generating at the point r of the 
detection device an interf erogram I (r, t) varying over 
time t, 

3 0 (f ) extracting a digital information item from 

the interf erogram I (r , t) , and 

(g) obtaining the coordinates of a point of 
measurement of the object to be imaged, to which the 
digital information item relates. 

35 

Thus the reference wave, which serves as local 
oscillator, need not pass through the object to be 
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imaged. This makes it possible to have a sufficient 
level of local oscillator, and to extract, with a 
better signal/noise ratio, useful information relating 
to the point of measurement for example for imaging 
5 purposes, in particular medical imaging. Moreover, this 
imaging method makes it possible to obtain a utilizable 
signal even with low acoustic or optical powers, that 
are for example compatible with the safety standards 
for imaged tissues associated with medical imaging. 

10 

In preferred embodiments of the invention, recourse may 
furthermore possibly be had to one and/or other of the 
following provisions: 

in the course of step (f ) , an acoustic 
15 component of the part of the scattered signal wave 
applied to the detection device is detected, this 
acoustic component being at a frequency corresponding 
to the sum of the frequency fi of the incident wave and 
of a harmonic of the acoustic frequency f A (fi ± H.f A , H 

2 0 nonzero integer) ; 

in the course of step (a) , said reference wave 
is generated at a frequency f R equal or substantially 
equal to the sum of the frequency f z of the incident 
wave and of said harmonic of the acoustic frequency 
25 fA- (£r « fi ± H.f A , H nonzero integer); 

in the course of step (b) , an acoustic wave is 
generated, focused at a focal point situated in the 
object to be imaged, and in the course of step (g) , the 
coordinates of the measurement point are obtained as 

3 0 being the coordinates of said focal point; 

steps (a) to (g) are repeated for various focal 
points of the acoustic wave in the object to be imaged, 
these various focal points being aligned along the 
first object direction; 
35 - in the course of a first iteration, steps (a) 

to (f) are performed for a first frequency f A of the 
acoustic wave and a first frequency f R of the reference 
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wave, in the course of at least a second iteration, 
steps (a) to (f) are repeated for a second frequency f' A 
of the acoustic wave and a second frequency f' R of the 
reference wave, these second frequencies being coded 
5 respectively with the first frequencies, the method 
furthermore comprising a step in the course of which: 

(f 1 ) at least one digital information item is 
obtained by decoding said digital information items 
obtained in the course of steps (f) of each iteration 
10 as a function of the frequencies used, 

and, in the course of step (g) , the coordinates of 
at least one point of measurement of the object to be 
imaged to which the digital information item obtained 
in the course of step (f 1 ) relates are obtained, by 
15 decoding said digital information items obtained in the 
course of steps (f) of each iteration as a function of 
the frequencies used; 

the sequence of following operations is 
performed: 

2 0 • a scan is performed of the frequency of the 

acoustic wave, which is focused on an interval of 
points with coordinates ( [U-Dx, U+Dx] , V, W) extending 
around the point with coordinates (U, V, W) along the 
first object direction, 
25 • a scan is performed jointly of the frequency 

f R of the reference wave in such a way as to keep f R 
substantially equal or equal to fi ± H.f A , H being a 
nonzero integer, 

• an interf erogram I (f A , V, W, r) associated 

3 0 with the set of points ( [U-Dx, U+Dx] , V, W) of the 

extended interval is recorded for each pixel r and for 
each f r equency f A ; 

• a ID frequency -» time Fourier transformation 
is performed, for each pixel r, according to the 

35 frequency f A of the interf erogram I (f A , V, W, r) , and 

© at least one interf erogram I (r) associated at 
least with a measurement point with coordinates (U' V, 
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W) is obtained by replacing the time obtained after the 
Fourier transform with the magnitude U 1 along the first 
object direction with the aid of the speed of 
propagation of the acoustic wave in the object to be 
5 imaged (U' possibly equalling U) ; 

at least steps (a) to (g) are repeated after 
having imposed a displacement of the vibration 
generating device relative to the object to be imaged 
along a direction not parallel to the first object 
10 direction of the object to be imaged; 

in the course of step (f ) , the complex 
amplitude E s ( r ) of the acoustic component is estimated 
on the basis of the interf erogram I (r , t) ; 

the detection device used is a monopixel 
15 detector, and, in the course of step (f ) , the digital 
information item is obtained as being the intensity of 
the field of complex amplitude E s ( r) scattered by the 
object ; 

the detection device used is a multipixel 

2 0 detector, and in the course of step (f ) , the digital 

information is extracted as being the sum over at least 
a part of the pixels r of the detection device of the 
intensity of the complex amplitude field E s (r) scattered 
by the object; 

25 - in the course of step (d) a spatial filtering 

device is used, in such a way as to limit, along at 
least one direction, the angular extent of the part of 
the scattered signal wave which is seen by each pixel 
of the detection device (it is thus possible to define 

3 0 a mean angular direction for the part of the scattered 

signal wave which is seen by each pixel of the 
detection device) ; 

a spatial filtering device comprising a 
diaphragm, of dimensions X along a first diaphragm 
35 direction and Y along a third diaphragm direction, and 
a lens of focal length L with object focus situated 
directly downstream of the object to be imaged is used 
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so as to limit the angular extent of the part of the 
scattered signal wave which is seen by each pixel of 
the detection device, and the reference wave applied to 
the detection device is globally a plane wave (the 
5 direction of application of the acoustic wave, of the 
incident wave, and the directions of the diaphragm are 
not necessarily linked) ; 

a spatial filtering device comprising a 
diaphragm of dimensions X along the first diaphragm 

10 direction and Y along the third diaphragm direction, 
disposed between the object to be imaged and the 
detection device at a distance L from the latter is 
used, so as to limit the angular extent of the part of 
the scattered signal wave which is seen by each pixel 

15 of the detection device, and the reference wave applied 
to the detection device is a spherical wave emanating 
from a source point situated in the plane of the 
diaphragm (the direction of application of the acoustic 
wave, of the incident wave, and the directions for the 

20 diaphragm are not necessarily linked) ; 

the reference wave and the scattered signal 
wave interfere on the detection device while forming a 
nonzero angle d Y , Oy being measured in the plane of 
incidence of these two waves on the detection device; 

25 - the detection device used is a multipixel 

detector, and the part of the acoustic component, of 
complex amplitude E s (r), which varies rapidly in space 
in the plane of the detection device is isolated (the 
detection device is in a plane quasi -orthogonal to the 

3 0 direction of the reference wave, and the part that one 
seeks to isolate corresponds to the components of the 
interf erogram I (r) , which vary rapidly in space and 
slowly in time) ; 

the detection device comprises pixels disposed 

3 5 as a matrix comprising rows along a first detector 
direction and columns along a third detector direction, 
and step (f) comprises the following steps: 
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(fl) for at least one row or column a ID Fourier 
transform is done along this row or column of the 
detection device to the space of wave vectors, of the 
complex amplitude E s (r ) , and a field TF X E s (k) , is thus 
5 obtained for this row or column, 

(f2) several zones of summation are defined in the 
space of wave vectors, 

(f3) the intensities of the field TF X E s (k) at each 
point k of at least one zone are summed in this zone, 
1 0 and 

(f4) the digital information item is extracted as 
being a linear combination of the sums thus obtained at 
each zone (this linear combination possibly comprising 
just one term) ; 

15 - the detection device comprises pixels disposed 

in a matrix comprising rows along a first detector 
direction and columns along a third detector direction, 
and step (f) comprises the following steps: 

(fl) a 2D Fourier transform is done of the complex 

20 amplitude E s (r ) , from the plain of the detection device 
to the space of wave vectors, and a field TF 2 E s (k) , is 
thus obtained, 

(f2) several zones of summation are defined in the 
space of wave vectors, 

25 (f3) the intensities of the field TF 2 E s (k) at each 

point k of at least one zone are summed in this zone, 
and 

(f4) the digital information item is extracted as 
being a linear combination of the sums thus obtained at 

3 0 each zone (this linear combination possibly comprising 
just one term) ; 

the angle 8 Y is about equal to 3Y/2L, in the 
course of step (f2) are defined a first zone of 
summation, the so-called central zone, a second zone of 

35 summation, the so-called left zone, and a third zone of 
summation the so-called right zone, and, in the course 
of step (f4) , the digital information item is extracted 
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as being a linear combination of the value of the sum 
of the left zone and of the sum of the right zone (this 
linear combination possibly comprising just one term) ; 
in the course of step (a) , 

• a laser source of wavelength X emits an 
emission wave, of frequency f L , 

• amplitude modulation means of the emission 
wave generate a carrier wave of incident 
frequency fi, and at least one amplitude 
modulation lateral band, which corresponds to 
a wave of frequency f R/ 

• a semiref lecting mirror, transmits a part of 
the lateral band wave and a part of the 
carrier wave forming the incident wave, and 
reflects a part of the carrier wave and a 
part of the lateral band wave forming the 
reference wave; 

in the course of step (a) , 

• a laser source of wavelength X emits an 
emission wave, of frequency f L/ 

• a first acousto-optical modulator 
transmits a part of the emission wave to 
form the incident wave on the object to be 
imaged, and moreover generates a first 
frequency shifted wave, the frequency of 
which is shifted by a value 6f lt possibly 
negative, with respect to the emission 
wave , arid 

• a second acousto-optical modulator 
intercepts the first frequency shifted 
wave and generates a second frequency 
shifted wave, the frequency of which is 
shifted by a value 6f 2 , possibly negative, 
with respect to the shifted wave, the 
second frequency shifted wave forming the 
reference wave, the frequency of which is 
thus shifted in frequency with respect to 
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the incident wave by a value 6f = 6fi+6f 2 , 
thus determining a known phase difference 
0i(t) between these two waves; (usually 6fi 
and 6f 2 are of opposite sign) ; 
in the course of step (a) , two independent 
laser sources, locked in phase by electronic slaving, 
generate the incident and reference waves, exhibiting a 
known phase difference <£i(t) between them; 
in the course of step (a) # 

• a laser source of wavelength X emits an 
emission wave, of frequency f L , 

• a semiref lecting mirror transmits a part of 
the emission wave to form the incident wave 
on the object to be imaged, and transmits a 
second part of the emission wave, 

• a first acousto-optical modulator intercepts 
the second part of the emission wave and 
generates a first frequency shifted wave, 
with frequency shifted by a value 6fi, 
possibly negative, with respect to the 
emission wave, and 

• a second acousto-optical modulator intercepts 
the first frequency shifted wave and 
generates a second frequency shifted wave, 
the frequency of which is shifted by a value 
6f 2 , possibly negative, with respect to the 
shifted wave, the second frequency shifted 
wave forming the reference wave, the 
frequency of which is thus shifted in 
frequency with respect to the incident wave 
by a value 6f=6fi+6f 2 , thus determining a known 
phase difference <f>±(t) between these two waves 
(usually 6f x and 6f 2 are of opposite sign) ; 

the object to be imaged is a biological tissue; 

the vibration generating device is used to 
obtain an acoustic cue of the zone of the object to be 
imaged, and the digital information item extracted in 
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step (f) is used jointly with said acoustic cue. 

According to another aspect, the invention relates to 
an installation for acousto-optical imaging of an 
5 object to be imaged (OBJ) comprising: 

means for generating an incident optical wave, 
of frequency fi, and a reference optical wave of 
frequency f R/ this reference wave being coherent with 
the incident wave and exhibiting therewith a known 
10 phase difference 0i(t), 

a vibration generating device for vibrating in 
a first object direction and at an acoustic frequency f A 
a zone of the object to be imaged, 

means for applying said incident wave to the 
15 object to be imaged, thus generating a scattered signal 
wave, 

a detection device, 

means for applying at least part of this 
scattered signal wave to the detection device, 

2 0 - means for applying the reference wave to the 

detection device without making it pass through the 
object to be imaged, thereby generating at point £ of 
the detection device an interf erogram I (r , t) varying 
over time t, and 

25 - means for extracting a digital information item 

and the coordinates of a point of measurement of the 
object to be imaged, to which the digital information 
item relates, from the interf erogram. 

30 In preferred embodiments of the invention, recourse may 
furthermore possibly be had to one and/or other of the 
following provisions: 

the installation furthermore comprises the 
following elements: 
35 • means for visualizing said digital 

information item relating to said point of 
measurement of the object to be imaged, and 
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• means for displacing the object to be imaged; 

the installation furthermore comprises a 
spatial filtering device situated downstream of the 
object to be imaged. 

5 

Other aspects, aims and advantages of the invention 
will become apparent on reading the description of 
several of its embodiments given by way of nonlimiting 
examples . 

10 

The invention will also be better understood with the 
aid of the drawings, in which: 

Figure 1 describes a first exemplary 
implementation of the method according to the present 
15 invention, 

- Figure 2 is a detailed diagram of an example of 
the device for generating two coherent waves according 
to the present invention, 

Figures 3 and 4 are graphs representing the 

2 0 phase difference between two coherent waves as a 

function of time, in a general manner and in a 
particular case, 

Figure 5 describes the first exemplary 
implementation of the method according to the present 
25 invention with the generating device of Figure 2, 

Figure 6 represents a mapping of the signal 
obtained, 

Figure 7 represents a second exemplary 
implementation of the method according to the present 

3 0 invention, 

- Figure 8 is a detailed diagram of another 
example of the device for generating two coherent waves 
according" to the present invention, 

Figure 9 represents a detail of a third 
3 5 embodiment of the invention, and 

Figure 10 represents the rear face of the 
diaphragm used in the third embodiment . 
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Figure 1 shows a device for generating waves GEN, which 
generates : 

an incident optical wave INC, of wavelength X, 
5 of frequency f lf applied to an object to be imaged OBJ, 
and 

a reference optical wave REF of frequency f R . 

The incident wave INC and reference wave REF are 
10 . mutually coherent and exhibit a known phase difference 
0i(t) . These optical waves may be emitted in the visible 
region, or possibly in the infrared or ultraviolet. 

The generating device GEN is adjusted in such a way 
15 that the reference wave REF is shifted in frequency 
with respect to the incident wave by a value equal to 
6f . 

There exist numerous possibilities for embodying such a 
20 generating device, and it is for example possible to 
use a generating device such as represented in 
Figure 2. This generating device comprises: 

a laser source LAS emitting an optical emission 
wave EMI at the frequency f L/ 
25 a first acousto-optical modulator MAOl 

transmitting a part of the emission wave EMI to form 
the incident wave INC at the frequency fi, and 
furthermore generating a wave DEC shifted by 5fi in 
frequency, and 

30 a second acousto-optical modulator MA02 

intercepting the frequency shifted wave DEC, of 
frequency f L + 8fi, so as to again shift this wave by a 
frequency 8f 2 to thus generate a reference wave REF at 
the frequency f R = f L + 8fi + 8f 2 , 8fi and/or 8f 2 possibly 

35 being negative (8fi and 8f 2 are usually of opposite 
sign) . 
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These acousto-optical modulators MAOl and MA02 are for 
example constituted by an acousto-optical cell of 
tellurium dioxide (Te0 2 ) , oriented along a given angle 
relative to the wave applied thereto, namely the 
5 emission wave EMI and the shifted wave DEC 
respectively, and vibrating under the action of a high- 
frequency generator, of frequency 8fi and 5f 2 
respectively, and transmit both a nondif f racted beam 
and a frequency shifted diffracted beam. The two 

10 acousto-optical modulators MAOl and MA02 performing the 
frequency shift of the reference wave REF with respect 
to the laser wave have for example frequencies of 
around 80 MHz, so that these two waves may be mutually 
shifted by a frequency shift 6f = Sfi + 5f 2 , which may 

15 be fixed between a few hertz and about a few tens of 
megahertz . 

As a variant, it is also possible to use another 
generating device (not represented) where the incident 
20 wave INC is extracted from the laser source LAS with 
the aid of a semitransparent mirror (or any other 
device) placed upstream of the first acousto-optical 
modulator MAOl, which then generates the shifted wave 
DEC on the basis of the part of the incident wave 

2 5 transmitted by the semitransparent mirror (or any other 

similar device placed between the laser source LAS and 
the first acousto-optical modulator) . 

It may also be possible to use two independent laser 

3 0 sources, locked in phase by electronic slaving, and 

generating incident INC and reference REF waves 
exhibiting a known phase difference (pi(t) between them. 

In these embodiments of a generating device, the laser 
3 5 source or sources LAS may for example comprise a laser 
diode stabilized by an external cavity with grating, 
emitting a laser wave EMI of wavelength X = 850 nm and 
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delivering a maximum power of 2 0 mW for a current of 
65 mA. It may be possible to use a power of said laser 
of about 15 mW, thus giving in this example, taking 
account of the optical losses, a power of about 7 mW of 
5 the incident wave INC on the object to be imaged OBJ. 

The two waves thus obtained, the incident wave INC at 
the frequency fi and the reference wave REF at the 
frequency f R/ are thus two coherent waves shifted in 
10 frequency by a value 8f = 8fi + 5f 2 . It is thus 
possible to control the phase difference q>i(t) between 
these two waves . 

Figure 5 represents the first embodiment of the 
15 invention with the generating device GEN of figure 2. 

The object to be imaged OBJ to which is applied the 
incident wave INC is a scattering object in respect of 
the optical waves, for example a sample of biological 
20 tissue. In the example considered, this sample may for 
example exhibit a thickness of about 2 0 mm in the 
direction of propagation of the incident wave INC. 

This sample may in particular be compressed between a 

2 5 front plate and a downstream plate, which are 

perpendicular to the direction of propagation of the 
incident wave INC, these two plates forming part of a 
sample holder (not represented) . 

3 0 The upstream plate is for example entirely transparent 

and made in particular of PMMA (Plexiglas ®) , while the 
downstream plate may for example be opaque and made in 
particular of black Bakelite. This* rear plate is 
pierced for example with a circular hole of diameter of 
3 5 about X = 20 mm. 

An acoustic coupling is effected between the obj ect to 
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be imaged OBJ and a vibration generating device TRANS 
situated outside the object. For this purpose, the 
object to be imaged OBJ and its sample holder may be 
installed at the center of a tank 1, for example 180 mm 
5 in diameter and 150 mm high. This tank may be furnished 
with plane glass windows 60 mm in diameter, 240 mm 
apart, and may possibly be filled with liquid serving 
to accomplish the acoustic coupling between the object 
to be imaged OBJ and the vibration generating device 
10 TRANS. 

If the object to be imaged OBJ may not be immersed, or 
for any other reason, other means, known to the person 
skilled in the art, may be used to guarantee the 
15 acoustic coupling. It is for example possible to mount 
the vibration generating device directly on the object 
to be imaged OBJ. 

The vibration generating device TRANS may be a 

2 0 transducer PZT possessing a frequency of acoustic 

vibration f A = 2.2 MHz, and a variable focal length 
fixed for example at 75 mm. It is excited for example 
by a sinusoidal signal produced at a maximum power of 
22 dBm on 50 Q, i.e. 6.3 Vcc, benign in respect of 
25 biological tissues. It may be possible to use a 
different maximum power, for example of the order of 
34 dBm or the like. 

This vibration generating device is oriented along a 
30 first direction of the object Xo, and emits an acoustic 
wave of frequency f A along this first direction of the 
object Xo. A zone (Dx, Dy, Dz) of the object to be 
imaged, centered on a point with coordinates (U, V, W) 
for which one seeks to obtain information is thus made 

3 5 to vibrate at the frequency f A by the vibration 

generating device TRANS. The extent (Dy, Dz) of the 
vibrating obj ect zone corresponds roughly to the 
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dimension of the focal zone of the acoustic wave, along 
the transverse directions y and z, that is to say in 
the plane normal to the direction of propagation Xo of 
the acoustic wave emitted by the transducer. In this 
5 plane, the focal zone is centered at (V, W) . U 
corresponds to the distance between the acoustic 
transducer and the focal point of the latter along the 
direction of propagation of the acoustic wave. Around 
this focal point, a certain zone, of variable extent Dx 

10 depending on the type of transducer used and the nature 
of the object to be imaged, among other things, 
vibrates at the acoustic frequency f A . Thus, the 
position and the orientation of the transducer and the 
position of its focal point determine a point of 

15 measurement of the object to be imaged OBJ, with 
coordinates (U, V, W) . 

The incident wave INC is applied to the object to be 
imaged OBJ, along a second object direction z Q/ possibly 
20 identical to the first object direction Xo, to form a 
scattered signal wave DIF which is scattered by the 
object in all directions. Inside the object to be 
imaged OBJ, a part of the wave passes through the zone 
(Dx, Dy, Dz) of the object to be imaged OBJ vibrating 

2 5 at the acoustic frequency f A . The movement of the points 

of the object that are likely to scatter generates a 
modulation of the phase of the scattered wave at the 
acoustic frequency f A . The vibration furthermore 
produces a modulation of the optical index of the 

3 0 medium (also at the frequency f A ) . These two effects 

give rise to the generation of an acoustic lateral band 
shifted in frequency by f A vis-a-vis the incident wave 
of frequency fi which passes through the medium. The 
scattered signal wave DIF therefore possesses a 
3 5 frequency shifted acoustic component of frequency 
foA = fi ± H.f A (where H = 1, 2, ... is the harmonic rank, 
in general H = 1) . The principle of acousto-optical 
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imaging then consists in detecting on a detection 
device DET this acoustic component of the signal wave 
scattered by the object, by making this acoustic 
component of the scattered signal wave interfere with a 
5 local oscillator of approximately similar frequency. 
This is performed according to the invention by using 
the reference wave REF not passing through the object 
to be imaged OBJ, as local oscillator for the detection 
device . 

10 

The detection device exhibits at least one detection 
cell in a plane x D , y D almost normal to the direction of 
observation z D (in the illustrative examples given in 
the figures, x D , y D/ z D correspond to Xq, y 0 , z Q ) . 

15 

It is possible to use a multipixel detection device 
DET, for example a digital CCD camera (12 bits) 
consisting of 1280 x 1024 square pixels of size 
d x = d y = 6 . 7 /-un. The CCD camera is chosen to have a 
20 sufficient quantum yield for a wave at 850 nm, for 
example 5% or more. The camera may be of the "full 
frame" type (non interlaced) or frame transfer type, 
with a detection frequency of for example f c = 12.5 Hz. 

25 A device is used to make the scattered signal wave DIF 
emanating from the object to be imaged OBJ interfere on 
the detection device DET with the reference wave REF. 
To perform effective detection on pixels of finite 
dimension it is necessary for these two waves to be 

30 almost colinear (forming an angle of 5° at most) . It is 
possible for example to use a semiref lecting pane or a 
beam splitter prism to guide the reference wave REF 
toward the detection device DET. 

3 5 We shall see that it may be useful to shift the 
direction of the reference wave REF angularly from the 
mean direction of the scattered signal wave DIF, by an 
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angle 0 Y . For th£ following description and in Figures 
1, 5, 7 and 9, the angular shift is performed along the 
direction y, the direction of the acoustic wave is x, 
the direction of observation is the same as that of the 
5 incident wave INC (direction z) and the matrix of the 
CCD detector is oriented along x and y, but these 
orientations are defined only by way of example. 

The scattered signal wave DIF and the reference wave 
10 REF are made to interfere on the detection device, and 
with the aid of the latter an interf erogram I (U, V, W, 
r, t) , is recorded, taken at the instant t, at the 
point r of the detection device, and which corresponds 
to the point of measurement with coordinates (U, V, W) 
15 of the object to be imaged, vibrating at the acoustic 
frequency f A . 

This installation may be coupled with a processing 
device CALC, capable of extracting from the time 

2 0 interf erogram recorded a digital information item 
relating to the measurement point, with coordinates (U, 
V, W) , this information item subsequently being 
displayable in an image of the object. This processing 
involves the calculation of the complex amplitude E S (U, 

25 V, W, r) of the acoustic component of the scattered 
field, frequency shifted by the acoustic vibration, on 
the detector. 

In a first nonlimiting embodiment of the determination 
30 of the complex amplitude E S (U, V, W, r) , the acoustic 
wave, of fixed frequency f A , is focused at the 
measurement point with coordinates (U, V, W) , and the 
frequency f R of the reference wave REF is chosen in such 
a way as to perform a N-phase demodulation, in 
35 particular a 4 -phase demodulation (N = 4) , as explained 
hereinbelow. 
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The measurement point, with coordinates (U, V, W) , of 
the object to be imaged OBJ vibrating at the acoustic 
frequency f A , the scattered signal wave DIF contains an 
acoustic component of frequency fi+f A . The processing 
5 unit CALC processes the interf erogram 1(11, V, W, r, t) 
which can vary over time by a four-phase demodulation 
in the following manner. The frequency shift 6f is 
firstly chosen between the incident wave INC and the 
reference wave REF such that 6f = f A +f c /4 = f A +3.125 Hz. 
10 The detection could however also be performed for any 
number N at least equal to 2 of phases, and the 
frequency of the reference wave adapted by taking 
6f = f A + f c /N. 

15 One then measures N interf erograms , each during a time 
Tint/ each interf erogram corresponding to a known 
distinct phase difference <f>± between the- incident wave 
INC and the reference wave REF. The reference wave 
being at the frequency f R = fi+f A +f c /N, and the acoustic 

2 0 component of the signal wave being at the frequency 
fi+f A , the phase difference between the reference wave 
REF and the acoustic component of the signal wave is, 
as represented in Figure 3 in a general manner and in 
Figure 4 for N equal 4, linear per interval as a 

25 function of time and goes from 0 to 211 over a time 
interval equal to N/f c - 

We thus detect N interf erograms Ii 1 (U, V, W, r, t) , 
I^fU, V, W, r, t) , In^U, V, W, r, t) corresponding 

30 to N distinct values <t>± of the phase difference, N being 
at least equal to 2, and in the case of Figure 4 being 
equal to 4 . 

This operation may be performed a number n at least 
35 equal to 1 of times, and for example for n=3, so as to 
obtain, for each known phase difference <f>±, n 
interf erograms Ii k (U, V, W, r, t) (i = l...N, k=l...n) . 
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For 12 interferograms, i.e. 3 cycles of 4 phases (N=4, 
n=3) , the total duration of the measurement is for 
example of the order of a second. 

5 

The 4 -phase demodulation calculation may be performed 
separately for each of the pixels r of the detection 
device. The processing device CALC carries out in fact 
for each pixel r of the detection device, the following 
10 operations: 

the integration of the intensity measured by 
the detector over the duration of a frame of the CCD 
camera (1/12.5 s) : Ii k (U, V, W, r) = (1/T int ) J I (U, V, 

W, r, t) dt taken between instants t = t ik and 
15 t = t ik + Tint- 

for each given phase <t>±, averaging of the n 
interferograms Ii k (U, V, W, r) (with k=l...n) detected for 
this phase to obtain N averaged interferograms Ii(U, V, 
W, r) . 

2 0 - a demodulation with N phases of the N 

interferograms measured so as to obtain the information 
regarding the complex amplitude E S (U, V, W, r) of the 
acoustic component of the signal wave DIF scattered by 
the object. For example, for n equal to 1, and in the 

25 particular case of Figure 4 where N is equal to 4 and 
where </>i = n/2, 0 2 = n, 0 3 = 3II/2, 0 4 = 211, the complex 
amplitude E S (U, V, W, r) of the acoustic component is 
proportional to (I 4 - I 2 ) + j di - I 3 ) (where j is the 
complex number such that j 2 =-l) , and the intensity 

30 associated with the object to be imaged OBJ also, since 
the reference wave is considered to exhibit a constant 
complex amplitude in space and in time. 

To determine the complex amplitude E s (U, V, W, r) of the 
35 acoustic component of the scattered signal wave DIF, a 
certain number of other techniques are commonly 
employed, known to the person skilled in the art. 
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A second mode of determination of the complex amplitude 
E s (U, V, W, r) of the acoustic component consists for 
example in using the so-called "frequency chirping" 
5 proceeding. In this case, instead of focusing at the 
magnitude U along the direction x an acoustic wave of 
fixed frequency f A/ it is possible to perform the 
sequence of following operations: 

a scan is performed of the frequency f A of the 
10 acoustic wave, which is focused on an interval of [U- 
Dx, U+Dx] extending around U; 

a scan is performed jointly with the first 
scan, of the frequency f R of the reference wave in such 
a way as to keep the condition f R = f z ± H.f A/ (where H 
15 is the harmonic rank, in general 1) (it will be noted 
that here the reference wave has the same frequency as 
the acoustic component that one seeks to detect) ; 

the interf erogram I(f A , V, W, r) associated with 
the set of points ( [U-Dx, U+Dx] , V, W) of the extended 
2 0 interval is recorded for each pixel with location r and 
for each frequency f A ; 

a ID frequency time Fourier transformation is 
performed, for each pixel with location r, according to 
the frequency f A of the interf erogram I(f A , V, W, r) , 

2 5 and 

the complex interf erograms I (U, V, W, r) 
associated with the various values of the magnitude U 
along the direction x are obtained by replacing the 
time obtained after the Fourier transform with the 

3 0 magnitude U along x with the aid of the speed of 

propagation of the acoustic wave in the object to be 
imaged . 

The reference wave being approximately a plane wave of 
35 constant amplitude, the complex interf erogram I (U, V, 
W, r) thus decoded is directly proportional to the 
complex amplitude E S (U, V, W, r) of the acoustic 
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component of the scattered signal wave DIF that one 
seeks to determine. 

There exist numerous variants of the "frequency 
5 chirping" technique which consist in replacing the 
coding along the x direction with a frequency coding 
according to f A and f R . These techniques, or others, may 
equally well be used to determine E S (U, V, W, r) in the 
context of this invention. 

10 

Provision may furthermore be made in this first 
embodiment for a spatial filtering device COL (Figures 
1 and 5) , which makes it possible to limit the angular 
extent of the part of the signal wave seen by each 
15 pixel of the detection device. We shall see that this 
device is useful for controlling the speckle grain 
dimension associated with the scattered signal wave 
DIF, in the plane of the detection device. 

20 To perform the detection it may be beneficial to make 
sure that the dimension of the speckle grains is 
adapted to the size of the pixels of the camera. This 
first condition corresponds to the so-called "anti- 
aliasing 11 condition . 

25 

Furthermore, the judicious choice of the geometry of 
this spatial filtering device will make it possible to 
isolate, according to the invention, the useful signal 
from the various components of noise. 

30 

This filtering device is for example constituted by a 
rectangular diaphragm 2, positioned for example 
directly downstream of the object to be imaged OBJ, 
perpendicularly to the direction of observation (and 
35 hence almost parallel to the detection device) , for 
example between the object to be imaged OBJ and the 
downstream plate of the sample holder, and elongated 
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along a direction. This diaphragm 2 may for example 
consist of two slender plates of aluminum sheet 0 . 5 mm 
thick separated by around Y=4 mm. The observable zone 
upstream is thus of quasi-rectangular shape of about 
5 X mm and Y mm along two directions perpendicular to the 
direction of observation. We may for example take z as 
direction of observation, and x and y as axes for the 
diaphragm. 

10 In the embodiment of the spatial filtering device COL 
represented in Figure 5, this device COL furthermore 
comprises a lens 3 placed between the tank 1 and the 
detection device DET. The object focus of the lens 3 is 
situated in the plane of the diaphragm (possibly taking 

15 account of the optical index of the liquid serving for 
the acoustic coupling of the vibrations generating 
device TRANS, if the scattered signal wave DIF passes 
through this liquid) . In the embodiment presented, a 
lens of focal length L=250 mm is used but other focal 

2 0 lengths may be suitable. 

This spatial filtering device reduces the angular 
extent of the part SIG of the scattered signal wave DIF 
which reaches the detection device DET, this possibly 
25 being useful for adapting the size of the speckle 
grains to the dimension of the pixels of the camera. 

Specifically, in the absence of any spatial filtering 
device, the signal wave DIF scattered by the object to 

3 0 be imaged OBJ may occupy a wide solid angle, of the 

order of n steradiants, and may be decomposed into a 
superposition of elementary plane waves with very 
different wave vectors Kg. Each wave vector Ks has, in 
the plane of the detector (x D , y D ) , two coordinates K x 
35 and K y . To obtain utilizable information, it is 
necessary for the density of the fringes corresponding 
to the spatial modulation of the signal of interference 
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of the scattered signal wave with the reference wave 
not to exceed the resolution of the matrix detector 
consisting of elementary detectors (so-called 
"anti-aliasing" condition) . In particular, for a plane 
5 reference wave REF, with wave vector Ko, with 
coordinates in the plane of the detector K x0 , K yo , this 
condition of density of the fringes gives rise to the 
"anti-aliasing" condition defined by: 

SKO = SINC (d x . (K x - K x0 )). SINC (d y . (K y - K y0 ) ) ~ 1 - 
10 Ea where the function SINC of a dummy variable xx is 
defined as being equal to 1 for xx = 0 and to sin xx/xx 
otherwise, where Ea is a decay factor quantifying the 
loss of contrast of the fringes that is related to the 
spatial integration of the detector, and where d x and d y 
15 represent respectively the characteristic dimensions of 
the elementary detectors of the detection device along 
the directions x D and y D . Thus, the measurement must be 
limited to an elementary angular field of the signal 
wave SIG, corresponding to a cone of angle (± Qf x ; ± a y ) 
2 0 around the direction of the wave vector Ko of the 
reference wave REF, the dimensions a x and a y of this 
elementary angular field having to be substantially 
less than or equal to X/2d x and X/2d y respectively so as 
to comply with said "anti-aliasing" condition, where X 

2 5 is the reference wavelength REF. 

The use of the diaphragm downstream of the object to be 
imaged OBJ, makes it possible among other things to 
eliminate the components of the signal wave not 
30 complying with this "anti-aliasing" condition. 

A judicious choice of the geometry of the spatial 
filtering device COL, and of the detection device makes 
it possible furthermore to isolate the useful signal 

3 5 from the various terms which appear in the signal 

resulting among other things from the analysis of the 
interf erograms I (U, V, W, r, t) . The matter is 
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discussed in the case of 4 -phase demodulation but a 
similar discussion could be had in the case of 
detection by "frequency chirping" , or other similar 
technique . 

5 

The interf erogram I (U, V, W, r, t) corresponds to the 
total intensity I T seen by the detection device, i.e. to 
the square of the modulus of the complex amplitude E 
(I T =|E| 2 =E.E* where E* is the complex conjugate of E) . 

10 To simplify the discussion we shall consider only a 
single acoustic component of frequency f A0 = fi + H.f A 
(with H = 1) . The complex amplitude E results from the 
sum of the amplitude E R of the reference wave REF, of 
the amplitude E T of the part of the scattered signal 

15 wave SIG at the frequency of the incident wave fi, and 
of the amplitude E s of the acoustic component of the 
signal wave of frequency f A o- We therefore have 
E=Ei+Er+Es - 

2 0 The total intensity I T/ which is used to extract the 

information item sought corresponds to the sum of 6 

terms (I = E.E* = (Ei+E R +E s ) . (Ei+E R +E s ) * ) : 

the term E x .Ei* corresponds to the interference 

between the part of the signal wave scattered at the 
25 frequency f x and itself, that is to say to the 

interference between the ordinary speckle and the 

ordinary speckle, 

the term E s -E s * corresponds to the interference 

between the acoustic component of the signal wave 
30 scattered at the acousto-optical frequency f A o = fi+fA 

and itself, that is to say to the interference between 

the acousto-optical speckle and the acousto-optical 

speckle, 

the term E R .E R * corresponds to the interference 
35 between the reference wave and itself, 

the term E R . Ei* corresponds to the interference 
between the reference wave and the ordinary speckle, 
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the term E^Es* corresponds to the interference 
between the ordinary speckle and the acousto-optical 
speckle, and 

the term E R .E S * corresponds to the interference 
5 between the reference wave and the acousto-optical 
speckle, which constitutes the term carrying the 
relevant information . 

The spatial filtering device COL makes it possible to 

10 reduce the angular extent of the wave emanating from 
the object which may behave at the level of the 
detection device as a quasi-plane wave. Such is the 
case in particular for the parts of the scattered 
signal wave having complex amplitude E x and E s . 

15 Additionally, the reference wave REF is, in the present 
embodiment, a plane wave. On account of the spatial 
filtering device the 3 diagonal terms E R .E R *, Ei.Ej* and 
E S .E S * (as well as the term Ei.E s *) vary slowly in space 
along the directions x and y of the detection device. 

20 Additionally the term Er.Ei* varies rapidly in time (at 
a frequency of around approximately f A ) and is averaged 
to zero on account of the low acquisition frequency of 
the detector. The relevant term E R .E S * is therefore easy 
to isolate by suitable digital processing (this making 

25 it possible to determine E s ) . If we choose a sufficient 
angular shift 0 Y between the reference wave and the 
scattered signal wave, this term is the only one to 
vary slowly over time, and rapidly in space along the Y 
direction. This Y direction corresponds to the 

3 0 direction of the width of the diaphragm, and to the 
direction Y D of the plane of the detector. 

One way of extracting the relevant information consists 
in performing a Fourier transformation of the complex 
35 amplitude E S (U, V, W, r) calculated above, along the x 
and y directions of the plane of the detector (or 
possibly the direction y only) . We then obtain a signal 
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TF E s (U, V, W, k) , k being the coordinate in the space 
of wave vectors. A mapping of the signal TF E s (U, V, W, 
k) obtained after Fourier transform is represented in 
Figure 6, which is an angular representation in the 
space of wave vectors. The various terms contributing 
to the interf erogram represented therein are discussed 
hereinafter. 

The N-phase demodulation (or the "frequency chirping") 
should in theory make it possible to completely 
eliminate the term for the interference between the 
reference wave and the reference wave (E R .E R *) , if the 
experiment were perfectly stable over time. This is 
never perfectly the case, and a fairly significant 
nuisance component therefore remains. However, this 
term varies slowly along the x and y directions of the 
plane of the detector, thus leading, in the space of 
wave vectors, on account of the bidimensional Fourier 
transform, to a narrow peak centered on the origin of 
coordinates (zone 4 of Figure 6) . 

The N-phase demodulation (or the "frequency chirping") 
should make it possible to eliminate the term for the 
interference between the ordinary speckle and the 
ordinary speckle (Ei.Ei*), if the experiment were 
perfectly stable over time, and if the speckle remained 
static and did not self decorrelate. This is never 
perfectly the case, and a fairly significant nuisance 
component therefore remains (zone 2 of Figure 6) . This 
term is even the dominant noise term for certain 
objects to be imaged in which the speckle does not 
remain static, (for example for certain biological 
tissues) . This term for the interference between the 
ordinary speckle and the ordinary speckle is, just like 
the term for the interference between the reference 
wave .and the reference wave, centered on the origin of 
the space of wave vectors. On account of the use of the 
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spatial filtering device defined previously, the 
amplitude of the ordinary speckle field has a finite 
angular extent, which corresponds to the interval 
[-Y/2L ; Y/2L] . As the term for the interference 
5 between the ordinary speckle and the ordinary speckle 
corresponds to the intensity of the field, that is to 
say to the product of the complex amplitude and the 
conjugate complex amplitude, it is necessary, in order 
to evaluate its extent in the space of wave vectors, to 

10 convolve the speckle field with itself in the space of 
wave vectors. This interference term therefore has an 
angular extent that is twice as wide as the field 
itself (interval [-Y 2 = -Y/L ; Y 2 =+Y/L] . Likewise, 
heightwise in the space of wave vectors, this term 

15 corresponds to the interval [-X 2 = -X/L ; X 2 =+X/L] , if 
the latter satisfies the aliasing condition in the 
optimum case (X/L = X/2d x ) . This noise thus exhibits a 
delimited pyramidal shaped envelope which is centered 
on the origin of the coordinates of the space of wave 

20 vectors . 

The interference between the acousto-optical speckle 
and the acousto-optical speckle (E S .E S *) is a second 
order diagonal term. Apart from its lower intensity, 
25 this term is no different from the term for the 
interference between the ordinary speckle and the 
ordinary speckle described above. 

The reference wave REF being at the frequency f R and the 
3 0 ordinary speckle being at the frequency f lt the term for 
the interference between the reference wave and the 
ordinary speckle (E R .Ei*) is at a frequency of about 
f R - f x about equal to f A i.e. about 2.2 MHz. This 
interference term is thus averaged to zero over the 
35 duration of acquisition of each image on account of the 
low acquisition frequency of the detection device and 
may therefore be neglected. 
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Likewise, the term for the interference between the 
ordinary speckle and the acousto-optical speckle 
(Ei.Es*)/ in addition to being a second order term, also 
5 possesses a frequency about equal to the frequency of 
the acoustic wave f A and may therefore be averaged to 
zero over the duration of acquisition of each image. It 
may therefore be neglected. 

10 The relevant term to be extracted from the interference 
between the signal wave SIG and the reference wave REF 
is therefore the term for the interference between the 
reference wave and the acousto-optical speckle (E R .E S *) . 
This term corresponds to zone 3 and to the interval 

15 [Y 3 _; Y 3+ ] of Figure 6. As this term corresponds to the 
interference of a speckle emanating from the diaphragm 
2 and of a plane wave, the interval [Y 3 _; Y 3+ ] has the 
same angular width as the angle of view of the 
diaphragm i.e. Y 3 . - Y 3+ = Y/L. This term evolves over 

2 0 time at the detection frequency f R -f A o (= f c /4 in the 

case of 4 -phase demodulation, and = 0 in the case of 
"frequency chirping") . It may therefore be measured by 
the electronic detection device. Additionally, as the 
reference wave is shifted angularly, the center of the 
25 zone 3 is shifted angularly by 0 Y (where the angle 0 y 
is, as defined previously, the angle between the 
reference wave REF and the signal wave SIG) . The choice 
of the angle 0 y defines the position of the center of 
zone 3, which is anyway of extent Y/L, when using a 

3 0 diaphragm of width Y. It is therefore appropriate to 

devise matters in such a way that zones 2 and 3 do not 
overlap so as to obtain, on the pixels of zone 3, 
useful information only. 

3 5 For example, we can choose 6 Y = 3Y/2L, this giving 
Y 3+ = 2Y/L, and Y 3 . - Y/L = Y 2 . Thus, the outer edge of 
zone 2 and the inner edge of zone 3 are superimposed 
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without the two zones overlapping. Neither is there 
then any blank part between zones 2 and 3 of Figure 6 . 
The signal of interference of the reference wave with 
the acousto-optical speckle has thus been separated in 
5 the space of wave vectors from the terms of the 
intensity of the reference and of the intensity of the 
speckle, which are situated at the center (zones 4 and 
2 of Figure 6) of the space of wave vectors. 

10 The extent of zone 3 being directly proportional to Y, 
there will be a temptation to increase the width of the 
diaphragm so long as the "aliasing" condition Y 3+ <X/2d Y 
is complied with. There is however a compromise to be 
made between on the one hand maximizing the area of 

15 zone 3 by increasing the width Y of the diaphragm, and 
on the other hand maximizing the integral of the 
product between area and effectiveness over this zone. 
This effectiveness corresponds to the loss of contrast 
of the interf erograms on account of the integration of 

2 0 the interf erograms over pixels of finite size. For 

adjoining pixels, an effectiveness factor is found 
according to a sine law similar to that introduced 
previously for the "aliasing". In this case, the 
effectiveness vanishes when the angle of "aliasing" 

25 X/2d Y is reached. For example, Y 3+ may be chosen to be 
appreciably less than the limit of "aliasing", for 
example Y 3+ = (2/3 ). X/2d Y approximately, so as to leave 
about a sixth of the area of the detection device 
inactive on each side, so that the effectiveness is 

30 sufficient. Y 3+ being equal to 2Y/L, this makes it 
possible to adapt, as a function of the laser used and 
of the detector used, the characteristics of the 
spatial filtering device. 

3 5 Figure 6 thus represents the mapping of the zones 

obtained in the space of wave vectors by the present 
invention after the bidimensional Fourier transform. 
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This map may be decomposed into a central column or 
zone 2, of extent [-Y/L, Y/L] , a left column or zone 1, 
and a right column or zone 3, Inside zone 2, zone 4 
represents the term for the interference between the 
5 wave of the reference and itself (E R .E R *) . To first 
order, if there is no acoustic signal, the noise 
observed on the image outside of zones 2 and 4 
corresponds to the "shot -noise" associated with the 
reference wave. In the presence of an acoustic signal, 
10 zone 3 represents the region of the space of wave 
vectors wherein is detected the useful signal according 
to the invention. According to the sign of the angle of 
incidence 6 Y , zone 3 could of course be situated on the 
left of Figure 6 . 

15 

The relevant information, which makes it possible to 
calculate the complex amplitude E s of the acoustic 
component, corresponds to the interference between the 
reference wave and the acousto-optical speckle (E R .E S *) . 

2 0 A digital information item relating to the measurement 
point (U, V, W) of the object to be imaged is for 
example obtained by summing the intensities calculated 
over the pixels of zone 3 ( | TF E s | 2 ) . The sum of the 
intensities of the pixels of zone 1 symmetric with zone 

25 3 with respect to the vertical axis of Figure 6 may 
serve as control. It is also possible to subtract from 
the sum of the intensities of the pixels of zone 3 the 
sum of the intensities of the pixels of zone 1, thereby 
making it possible to reduce any systematic measurement 

30 errors. The sum of the intensities measured over zone 2 
also carries information, since it characterizes the 
level of noise associated with the decorrelation of the 
speckle. 

35 It is also possible to use conventional means of image 
processing, to average, after 4-phase demodulation, 
Fourier transform and calculation of the intensity, the 
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intensities of the pixels for example 8 by 8, 16 by 16, 
or 32 by 32, so as to display a control mapping of the 
signal obtained for the measurement point, with 
coordinates (U, V, W) (to verify for example the 
5 positioning of the various zones) . 

We may also only carry out for each row of the 
detection device, a unidimensional Fourier transform 
(along y for the examples considered) , in which case 
10 the signal TF X E s (U, V, W, k) is not exactly the one 
represented in Figure 6, but zone 4 extends over the 
entire X-Y extent of this figure. The remainder of the 
processing remains valid. 

15 We thus obtain the sought-after digital information 
item, relating to the measurement point, with 
coordinates (U, V, W) of the object to be imaged OBJ, 
in the case of 4 -phase demodulation. For "frequency 
chirping", or the like, the digital information item 

20 thus obtained relates to the acoustic frequency f A used, 
and we obtain the coordinates and the digital 
information relating to various points of the interval 
[U-Dx; U+Dx] by the means described above. 

25 The choice of the angle 6 Y makes it possible to properly 
separate the signals obtained, but the appropriate 
positioning of the device making it possible to use a 
given angle 0 Y may require a control step. In this 
control step we obtain an image, for example by 

30 suppressing the acoustic wave and adjusting the 
frequency f R of the reference wave, so as to detect the 
component of the scattered field at the frequency f x . We 
choose for example f R « f z + f c /4 and we perform a 
4 -phase demodulation. The lens is positioned precisely 

3 5 in such a way that a sharp image of the zone seen 
through the diaphragm 4 is obtained by Fourier 
transform of the signal detected in the plane of the 
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detection device. 

As seen previously, in a Fourier transform calculation, 
the edge of the calculation matrix corresponds to the 
"aliasing" condition. The control step makes it 
possible to ensure, for example, that the outer edge of 
the zone containing the useful signal is not too close 
to the edge of the calculation matrix, and/or that the 
zone containing noise and the zone containing the 
useful signal are in contact at the level of the inner 
edge of this zone, but do not overlap. 

Figure 7 represents a second exemplary implementation 
of the method according to the invention in which the 
spatial filtering device is not used. Specifically, for 
low laser intensities, and if the speckle decorrelation 
term is not too large, the "shot -noise" becomes the 
dominant noise. Such is the case when with the first 
exemplary implementation (Figure 5) , the values of the 
sums, measured over zones 1 and 2, are about the same. 
This configuration is also useful for other 
applications where the acoustic wave is of sufficient 
power to reach a sufficient acousto-optical conversion 
yield. The signal associated with the acoustic 
component may then be larger than the speckle 
decorrelation noise. 

It is possible to dispense with the spatial filtering 
device since, even if the term for the interference 
between the ordinary speckle and the ordinary speckle 
(Ei.Ei*) can then no longer be separated spatially from 
the relevant term (E R .E S *) / it remains very much less 
than the "shot-noise" of the reference wave. It is also 
possible in this case to choose a zero angle d Y , for 
example with the aid of a semiref lecting pane. Only the 
term for the interference between the reference wave 
and itself (E R .E R *), which is highly centered on the few 
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pixels of zone 4 of Figure 6, constitutes a noise term 
to be filtered. This filtering is carried out simply by 
digitally eliminating the components of the space of 
wave vectors round about (0, 0), excluding zone 4 from 
5 the summation over the pixels of the intensities of the 
signal. It may even be that this noise term is 
negligible, and it is then possible to compute the sum 
of the intensities over the pixels in the real space 
rather than in the space of wave vectors (and it is no 
10 longer necessary in this case to perform a Fourier 
transformation) . 

Figure 8 represents another generating device used to 
implement the method according to the invention in the 

15 particular case of the second example of Figure 7. The 
emission wave EMI, of wavelength X, emitted by the 
laser LAS, is intercepted by an amplitude modulation 
device MA. This device generates a carrier wave POR, of 
frequency f lt and two amplitude-modulated lateral bands 

2 0 LATMOD and LATMOD 1 , of frequency f R . These three waves 
are applied to a semi transparent pane, which transmits 
a part of each of these waves applied to the object to 
be imaged OBJ. The transmitted part of the carrier wave 
constitutes the incident wave INC. Furthermore, the 

2 5 semitransparent pane reflects a part of each of these 

waves towards the detection device, the reflected part 
of LATMOD constituting the reference wave REF. 

This device is only adapted to the second embodiment, 

3 0 where no spatial filtering device is used, since in the 

case of an intensity modulation lateral band, the 
spatial filtering does not make it possible to separate 
the various interference terms. 

3 5 Figure 9 describes a detail of a third embodiment of 
the present invention and of the spatial filtering 
device. Numerous characteristics have already been 
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described when presenting the first embodiment, and 
will therefore not be described again. In this third 
embodiment, the reference wave REF is no longer a plane 
wave, but a spherical wave of frequency f R . Such a 
5 spherical reference wave may for example be obtained, 
on the basis of the plane reference wave generated by a 
generating device GEN described previously, by focusing 
the plane reference wave, with the aid of a lens 5, 
onto a small mirror 4 situated for example in the plane 

10 of the diaphragm 2. The reference wave REF may moreover 
arrive on the mirror 4 at an angle 0 Y so that it is 
reflected and arrives at the detection device DET while 
forming an angle d Y with the signal wave SIG. In this 
embodiment, one does not necessarily use a convergent 

15 lens 3 to collimate the signal wave heading for the 
detector DET, specifically, the divergent nature of the 
reference wave fulfils a similar role. 

Furthermore, the diaphragm 2 can be taken out of the 
20 sample holder of the object to be imaged OBJ, and can 
be placed between the tank and the detection device 
DET, or possibly can be fixed onto the front face of 
the tank 1. This provision may be beneficial if the 
acoustic coupling between the transducer TRANS, of 
25 acoustic frequency f A , and the object to be imaged OBJ 
is effected by a tank filled with water 1, and when the 
latter is situated on the light path travelled by the 
reference wave REF before reaching the mirror 4. (This 
provision is also valid for the other embodiments where 
3 0 a diaphragm is used) . 

The processing of the interf erograms measured is then 
identical to the processing setforth previously, where 
L now represents the distance between the diaphragm 2 
3 5 and the detection device DET. 

Figure 10 represents the diaphragm 2 furnished with the 
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mirror 4, visualized from the detection device DET. The 
diaphragm 2 exhibits a slit, for example central, of 
width Y and of height X, To obtain an image such as 
represented in Figure 6, a mirror 4, for example 
5 circular, is placed midway up the height of the slit. 
The mirror 4 is moreover shifted laterally with respect 
to the slit of the diaphragm 2. This lateral shift is 
related to the shift of zones 2 and 3 of Figure 6 that 
one seeks to effect, once the calculations have been 

10 performed by the processing device CALC. By choosing 
for example an angular lateral offset equal to 6 Y 
between the source point of the reflected reference 
wave and the edge closest to the mirror of the slit of 
the diaphragm 2, seen from the detection device, we 

15 obtain zones 2 and 3 of Figure 6 juxtaposed but not 
superimposed. In this case, care will of course be 
taken to choose a mirror 4 of dimensions such that the 
mirror does not encroach on the slit of the diaphragm 
2. Thus, it is easy to separate zones 2 and 3 of Figure 

20 6. 

This embodiment makes it possible moreover to easily 
find an optimal position of the detector DET making it 
possible to obtain information of good quality for the 

2 5 point studied (U, V, W) of the object to be imaged OBJ. 
Specifically, by increasing L, that is to say by 
distancing the detection device DET from the diaphragm 
2, the angle between the signal wave SIG and the 
reference wave REF decreases, so that the "aliasing" 

30 condition is better complied with. On the other hand,, 
the intensity of the signal decreases with this 
distance L. The compromise to be made between the 
effectiveness of detection and the size of the zone 3 
of Figure 6, already mentioned, is effected simply by 

35 the positioning of the detection device DET with 
respect to the diaphragm 2 . 



f 



WO 2004/085978 - 37 - PCT/FR2004/000640 

In this embodiment, the optimal positioning of the 
detector is carried out independently of the relative 
positioning of zones 2 and 3 of Figure 6, which, 
itself, is carried out by the appropriate placing of 
5 the point of focusing of the reference wave (center of 
the mirror 4 on the diaphragm 2)\ 

By implementing the method described previously, 
according to any one of its embodiments represented in 

10 Figures 1 to 10, we have obtained digital information 
relating to the measurement point, with coordinates (U, 
V, W) of the object (in the case of N-phase 
demodulation) or relating to each point of the interval 
( [U-Dx,U+Dx] , V, W) of the points of the object (in the 

15 case of "frequency chirping" or other similar 
procedure) , that is to say pointwise or ID information. 

To obtain an image of the 2D or 3D object, it is 
required to repeat the whole set of operations 

20 described above after having displaced the position of 
the focal point of the acoustic wave along the other 
directions (x or y or z for the N-phase demodulation) 
(y or z for the "frequency chirping"). The displacement 
of the position of the focal point may be achieved, 

25 either by displacing the transducer (along the 
direction x, y or z considered) while preserving the 
acoustic coupling between the transducer and the object 
to be imaged, or by using several transducers with 
different focal points, for example. 

30 

It may be beneficial to displace the spatial filtering 
device COL (if one is being used) , the detector DET, or 
even the device used to make the reference wave 
interfere on the detector, and all or part of the 
3 5 generating device GEN conjointly with the displacement 
of the position of the focal point i.e. of the point of 
measurement with coordinates (U, V, W) . Alternatively, 
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it may possibly be preferable to displace only the 
object to be imaged in the opposite direction so as not 
to have to perform a similar movement on the whole of 
the remainder of the installation. 

5 

As an acoustic transducer is used, any device according 
to the invention may furthermore be coupled to a 
conventional acoustic imaging device. The coupling of 
these two devices thus makes it possible to obtain from 
10 the object to be imaged OBJ purely acoustic information 
provided by the acoustic transducer by a conventional 
technique of ultrasound imaging, and optical 
information provided by the device according to the 
invention. 



